Two hadron production in e + e annihilation to 
next-to-leading order accuracy* 



D. de Florian, L. Vanni 

Departamento de Fisica, FCEYN, Universidad de Buenos Aires, 
(1428) Pabellon 1 Ciudad Universitaria, Capital Federal, Argentina 



Abstract 

We discuss the production of two hadrons in e + e~ annihilation within the framework 
of perturbative QCD. The cross section for this process is calculated to next-to- 
leading order accuracy with a selection of variables that allows the consideration of 
events where the two hadrons are detected in the same jet. In this configuration we 
contemplate the possibility that the hadrons come from a double fragmentation of a 
single parton. The double-fragmentation functions required to describe the transition 
of a parton to two hadrons are also necessary to completely factorize all collinear 
singularities. We explicitly show that factorization applies to order a s in the case of 
two-hadron production. 
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1 Introduction 



The production of one hadron in e + e~ annihilation has been studied in much detail in perturbative 
QCD [l]. The corresponding cross section for the process e + e~ — > J*(Q) — >• H(P) + X is usually 
expressed as a function of the variable 

representing the energy fraction carried by the hadron. In this case the cross section can be 
written as a convolution of the (perturbative computable) partonic cross section a 1 and the (non- 
perturbative) fragmentation functions Df(x) giving the probability of finding a hadron in the 
parton with momentum fraction 

£ = E**«tf- (2) 



The cross section has been computed to next-to-leading order (NLO) accuracy in pQ and 
to next-to-next-to-leading order (NNLO) accuracy in [2]. Furthermore, several analyses of the 
available data have been performed in the last years and, as a result, fragmentation functions for 
several hadrons have been extracted with very good precision. 

Higher order QCD corrections (NLO) to the cross section for the production of two hadrons 
in e + e~ annihilation have been computed in pQ in the particular case when the two hadrons H% 
and H 2 are selected form different parton jets. While a symmetric extension of the one-hadron 
case to two hadrons would correspond to expressing the differential cross section in terms of the 
momentum fractions of each hadron defined by 

_ 2Pi ■ Q _2P 2 -Q 
the authors in pQ introduced a different set of variables 



While z in Eq.Q coincides with z±, the momentum fraction of hadron Hi, the second variable 
u depends on both the momentum fraction of hadron H 2 and the angle 9 12 between the hadrons 
observed from the center of mass system as 

u = z 2 ~(1 - cos 6>i 2 ) , (5) 

such that u is approximately zero when the angle between the hadrons is small. Therefore, 
configurations where both hadrons are in the same parton jet corresponds to u w 0. Consequently, 
by considering events were z and u are not too small one can ensure that the two hadrons are 
produced from the hadronization of different partons and the cross section can be reduced to 
the product of the fragmentation functions Df associated to each hadron pQ. In this way, the 
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possibility of a double fragmentation from a single parton is excluded and the expression for the 
cross section gets simplified. 

In this work we are interested in extending the calculation for the two-hadron cross section 
in the full phase space, including the configurations were both hadrons are produced collinearly. 
In order to be able to consider those events we will express the cross section in terms of the 
momentum fractions in Eq. (J3J) • 

With the use of these variables it is possible to contemplate simultaneously two extreme con- 
figurations: the first one corresponds to the case when the two hadrons are produced in opposite 
directions (or at least with a clear angular separation) and therefore belonging to different jets 
(Figure^,). Hadrons in this configuration can only be originated from the fragmentation of dif- 
ferent partons. The second one corresponds to the case of both hadrons produced in the same 
direction, such that they are detected in the same jet (Figure UJp). In the last case, hadrons could 
be originated from the fragmentation of two collinear partons or by the double fragmentation of 
the single parton. The price one has to pay to fully account for the second configuration is the 
introduction of a new set of non-perturbative phenomenological functions [3J HJ |5] to describe the 
possibility of the transition from a single parton to two hadrons. 




H 2 

(a) (b) 



Figure 1: (a) represents hadrons in the first kinematical configuration, (b) represents hadrons in 
the second configuration, belonging to the same jet. 

One needs to introduce, then, the double-fragmentation functions DDp lH2 (xi,X2) 2 , as the 
probability that a parton p fragments into the hadrons H\ and H 2 with energy fractions x\ and 
x 2 . The cross section for the production of two hadrons in e + e~ annihilation can therefore be 
written in the following schematic way 

d -^- = Y 0* ® D?> ® 2?f a + W ® DD?*>, (6) 

dz x dz 2 ^ 3 ^ 

ij 1 

where a 13 is the partonic cross section for the production of partons % and j and a 1 the cross 
section for parton i. The cross section is separated in two terms corresponding to the contribution 
of the mechanisms responsible for two-hadron production: single fragmentation of two partons, 
and double fragmentation of a single parton. 

At leading order, the first term only contributes to the first configuration, since the two partons 
that undergo hadronization are produced back-to-back. At next to leading order there is one extra 

2 We have slightly modified the original notation introduced in [21 El f° r the double-fragmentation functions to 
make more noticeable the difference with the usual ones 
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parton which could be emitted collinearly to one of the others, giving also origin to hadrons in 
the second configuration. Therefore, at order a s and beyond, hadrons in the second configuration 
could be originated from any of the two fragmentation mechanisms, being not possible to separate 
the contribution of each term in Eq.(JHJ), unless an additional (unphysical) scale is introduced. 
Only the sum of both contributions has physical sense. 

The presence of collinear partons at order a s gives origin to collinear singularities in the cross 
section, which are manifested in the form of poles in e = (4— N)/2 when dimensional regularization 
is used. By means of the usual redefinition of the the Df functions, singularities due to collinear 
partons that give origin to hadrons in the first configuration can be absorbed. However, there 
appear singularities corresponding to hadrons belonging to the second configuration, originated 
from collinear partons emitted in the same direction. Since at lowest order the Df functions 
only participate in processes associated with the first configuration, such singularities can not 
be absorbed in the single fragmentation term. We will show that with the redefinition of the 
DDf lIi2 functions in the double fragmentation term all singularities are factorized. In this sense, 
the role played by the DDf lH2 functions in e + e~ annihilation is similar to the one of fracture 
functions in DIS processes |H1 El [HE [TTj . For a formal point of view, it is possible to interpret 
the double-fragmentation functions as the time-like version of fracture functions. 

Double-fragmentation functions DDf lH2 fulfill sum rules in analogy to the sum rules for the 
usual fragmentation functions pQ. Momentum conservation requires 



#1 



(] a HiH 2 j H 2 



being Q the initial total momentum, where the right hand side is proportional to the total free 
momentum available for the production of hadron H 2 . In particular, energy conservation implies 

/ d Zl Zl DDf lH »(z u z 2 ) = (1 - z 2 )D^(z 2 ), (8) 
relating the second moment of the double-fragmentation function to the single one. 



2 Two-hadron production in e + e 

In order to formalize the convolution products in Eq.© we define the partonic energy fractions 
associated to each fragmentation mechanism. For the single fragmentation term two partons 
fragment independently with a fraction of the parent parton energy given by 

* = ^, (»> 

with pi the momentum of the parton % = 1,2. At leading order both variables are fixed to one 
since no extra gluon radiation is allowed. 

The convolution product in the single fragmentation term of Eq.© is expressed in terms of 
a double integral in x\ and x 2 with integration intervals determined by the kinematical region 
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allowed for the partonic process. This implies 



< X\ < 1 and Z\ < Xi 
1 — x\ < X2 < 1 and 22 < £2 



(10) 



The integration zone for the single fragmentation term has to be divided into the two regions 
A and B indicated in FigEl 




Figure 2: Integration regions for the variables x\ and x 2 in the single fragmentation term. 



In the case of the double fragmentation term only one parton fragments. We define the partonic 
variable as usual by 

2p-Q 



x 



Q 2 
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with p being the momenta of the fragmenting parton. With this, it is possible to write the second 
term of Eq.fjHJ) as a single convolution product with integration limits coming from the request 
zi/x + z 2 /x < 1. 

Using those conditions we can express Eq.fJHJ) as 



da H ^ 
dz\dz 2 



1-Z2 



dxi f 1 dx2 do % i 



+ C 1 dxi C 1 dx 2 dcr lj dHi /Zi\ h s 

„Ji-xa Xi J Z2 x 2 dx 1 dx 2B 1 Vxi/ 3 v 



£2 

X 2 



1 dx da 1 

Z1+Z2 



^ I x 2 dx % 



X X 



(12) 
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where to NLO accuracy 



da 1 da i{0) a s da iW 
dx dx 2tt dx 

do*' do ij{0) a, do ij{1) 



dx\dx<iK dx\dx 2K ~^ ^ n dx\dx%^ ^ 

and K = A,B indicating the integration zone in the single fragmentation term. In Eq.fjl2j) we 
have considered the case when 1 — z 2 > z±. This implies z\ + z 2 < 1, which corresponds to the 
kinematical region where the double fragmentation mechanism can also contribute. If z\ + Z2 > 1 
the cross section is reduced only to the second term of Eq. (JT2J) . 

Some of the partonic cross sections obey symmetry relations that allow to reduce the number 
of independent quantities to be computed. Due to invariance under charge conjugation 

do iq do iq da qi do qi , . 

(14) 



dxidx 2 dxidx 2 dx\dx 2 dx\dx 2 

To NLO accuracy it is necessary to obtain only three different partonic cross section do qq j dx\dx 2l 
do q9 jdx\dx 2 and do 9q / dx\dx 2 . 

At leading order the only non- vanishing terms are 3 

do*<® 



e 2 o S(l - x) 



dx 

do qqi0) 
dx\dx 2K 



e 2 a o 6(l-x 1 )5(l-x 2 ). (15) 



1 Anai 

where cr = -^i 



The partonic cross section at order a s is obtained by evaluating the real and virtual diagrams 
indicated in Fig|3] and integrating over the phase space of the final partons expressed in terms of 
xi,x 2 , such that do^ = doR + day. We compute the metric and longitudinal contributions to the 
partonic cross section obtained, as usual, by replacing the sum over the polarization states of the 
virtual photon by the corresponding projectors 

= -9,,, (16) 

p - = j£w p2 ^ (17) 

The longitudinal contribution has been calculated projecting in the direction of hadron H 2 . 
In the following, we will present in detail the results for the metric contribution, since at NLO 



3 In this work we restrict the analysis to the case of pure 7* exchange. The extension to Z-boson production 
can be easily performed 
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Figure 3: Virtual and real diagrams contributing to order a s 



singularities of interest occur only on that projection of the cross section. Using dimensional 
regularization ^21 EH! we obtain for the real part 



da 



(M) 
R 



2 a s n 

e q a 07r: c F 











r(2-e) 





with 



F(x 1 ,x 2 ) 



'1 -e 



rf> 2 I ry^ 2 

2 X 1 -|- X 2 



2e(l-e) 



x x 2e x 2 2e F(xi, x 2 )dxidx2, (18) 
(2 — 2xi — 2x 2 + ^1^2) 



(1 - Xi)(l - x 2 ) ' ' (l-x 1 )(l-x 2 ) 
and /i being the dimensional regularization scale. The result for the virtual contribution is 

4Wr(i + e)r 2 (i -e) 



(19) 



day = e q a -C F{ Q2 



r(i-2f 



+ 71' 



5(1 — xi)S(l — x 2 )dxidx 2 .(20) 



In the previous equations we have labeled the quark as parton 1 and the anti-quark as parton 
2. The virtual part does not exhibit singularities beyond those already regularized in the form of 
poles in e. In the real part the divergences appear when the denominators of the function F(xi, x 2 ) 
vanish. These infrared divergences can be regularized by means of the usual + prescription, which 
can be easily implement by multiplying and dividing by (1 — Xi) 1+e and considering the distribution 



, 1 



-1-e 



i 6{1 ~ Xi) + ah-) 



log(l - Xj 

1 X/, 



+ 0{e 2 ) 



(21) 



-[0,1] 



where F(x) + [ at b] is defined as usual by 

dxf(x) F(x) + \ a> 



dx[f{x)-f{b)]F{x) 



The range of integration is indicated as a subscript; furthermore, the subtraction point is 
underlined. 

For do~ q ^ 1 ' dx\dx 2 the singularities of the function F(x\, x 2 ) occur at x 2 — 1 in zone A, and at 
X\ — 1 and x 2 = 1 in zone B. Applying the + prescription as indicated, the following expression 
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is reached 



da qq 



(i) M 



dX\dx 2K 



P qq ( Xl ) log ( 5(1 - x 2 ) + P qq (x 2 ) log ) 5(1 - Xli 



+ -J> qq {x x )8{\ - x 2 ) + jP qq (x 2 )5(l - Xl ) + fM{x u x 2 ) 



(22) 



where 1/e = -l/e(47r)T(l - e)/r(l - 2e) = (-1/e + 7s - log4vr) + 0(e), i%-(x) are the usual 
Altarelli-Parisi splitting kernels ^3] and the functions ftf K are presented in the Appendix. 

The qg partonic cross section da qg ^ / ' dx\dx 2 can be obtained from the qq one relabeling the 
parton indexes by a x 2 —> 2 — x\—x 2 substitution in the matrix elements. In this case F(x x , x 2 ) — > 
F(xi, 2 — xi — x 2 ) which develops singularities at x± = 1 in zone A, and at x\ + x 2 = 1 in zone 5. 
Proceeding like in the previous case and ignoring terms containing distributions without support 
in the analyzed zones, we obtain 



da qa{l) 
dx\dx 2K 



(M) 



PZM log ( S( Xl +x 2 -l)+ P* q (x 2 ) log ) 5(1 - Xl ) + 



+ -P q gq (x l )8(x 1 + x 2 - I) 



jP^2)S(l-X 1 )+f^(x 1 ,X 2 ) 



(23) 



where the functions P^ are the real LO Altarelli-Parisi kernels, with the index k labeling the third 
particle in the vertex % — >■ jk. 

In the gq partonic cross section da 9q ^ /dx x dx 2 , X\ is assigned to the gluon and x 2 to the quark. 
Performing the substitution x 2 <-* X \ in the matrix element used in da q9 ^ /dx x dx 2 , F(x±, 2 — x x — 
x 2 ) — > F(x 2 , 2 — xi — x 2 ) which develops singularities at x 2 — 1 and x x + x 2 = 1 in zone A, and at 
x 2 = 1 in zone B. The singularities in zone A will give origin to two different distributions, one 
associated to the singularity at x 2 = 1 and another to the singularity at x\ — 1 — x 2 . The result 
can be expressed as 



da 9q{l) 
dx\dx 2K 



(M) 



P q q ( Xl ) log ( 5(1 - x 2 ) + P 9 qq (x 2 ) log (^) 8( Xl +X2-1) 



+ \P q gq (x X )5{l - X 2 ) + jPy&Wxi + X 2 - 1) + /^(X!, x 2 ) 



(24) 



This result completes the presentation of the partonic cross sections that participate in the 
single fragmentation term. 

For the double fragmentation term, the cross sections for the production of a single parton 
do 1 jdx are required. These are exactly the same as the ones appearing in one-hadron production 
PP. As a cross-check of our calculation we have re-obtained those coefficients by applying the 
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momentum conservation relation in Eq.(J7J) resulting in 



doc doc 
do*® {M) 



P qq {x) lo£ 



/i 2 



(AO. 



X 



dx 



2e 2 a 



P gq {x)\og 



X) 



with f^ M ) given in the Appendix. 



(25) 



As indicated above, the longitudinal part does not contribute to the singular structure of the 
cross section to NLO accuracy. The corresponding NLO corrections to the single ftp and double 

f^ L ) fragmentation mechanism are shown in the Appendix. 



3 Factorized fragmentation functions 



The factorization of the bare fragmentation functions Df at NLO is done in the MS factorization 
scheme in the standard way PQ. The expression for the bare functions in terms of the factorized 
ones at the scale M 2 is 



D?{z) 



1 du 
u 



5(1 - u)8ij + 



2^ 



log 



M 2 



PiiM 



Df {NLO) (-,M 2 n 



. u 



(26) 



where the factorized distributions are labeled by the upper index (NLO). 



It is easy to notice that not all the singularities in the partonic cross-section are canceled 
after the factorization of the fragmentation functions is performed. Singularities belonging to 
do q9{1) /dx 1 dx 2 and do 9q(l) / 'dx 1 dx 2 in zone A still remain. They correspond to terms with 1/e 
poles proportional to 5(x\ + x 2 — 1), arising from the hadronization of a gluon being emitted 
collinear to a quark (or anti-quark) that also undergoes hadronization, giving as a final product 
two hadrons in the same jet. Those singularities clearly cannot be absorbed by the factorization 
of the Df functions, since a configuration with two collinear hadrons is not allowed in the single 
fragmentation term at the lowest order. However, this is exactly the configuration corresponding 
to the double fragmentation term, indicating that these singularities could be absorbed by the 
appropriate factorization of the DD^ lH2 functions. 

The expression for the bare double-fragmentation functions DDj IlH2 (x, y) in terms of the NLO 
factorized ones DDf lH2<yNLO \x ) y ) M 2 ) can be obtained by requiring that all remaining collinear 
singularities in the partonic cross section are absorbed into the factorized distributions. The 
expression in the MS factorization scheme, valid to 0(a s ), is 



DD 



HiH 2 . 



x,y) 



x+y 



du 
u 2 



5(1 — u)5ij + 



2tt 



log 



JUL 

M 2 



P jt (u) 



DDf^ NL °)(*,y,M 2 
J \u u 



The factorization relation in Eq. (|27|) contains two terms with different origins. The first one just 
relates the factorized and bare double-fragmentation functions, and corresponds to the standard 
factorization procedure for the emission of collinear partons in the double fragmentation part of 
the cross section, exactly as it occurs for one-hadron production. The second 'inhomogeneous' 
term relates the single and double-fragmentation functions and is needed to absorb the remaining 
singularities discussed above. 

Rewriting the bare distributions in terms of the factorized ones in Eq. (jl2j) . and fixing M 2 = Q 2 , 
we obtain the final NLO expression for the factorized cross section for the production of two 
hadrons as 



da H ^ 
dz\dz2 



f 1 Z2 dxi f 1 dx-2 v-^ 2 J "■- ri u.. 



X\ Ji- Xl x 2 



D H l{ NLO) f Q2) D ^(NLO) f *2 q2 \ + ^NLO) f H Q 2) MNLO) f * q2 

q \xi J \x 2 / Vx 2 / \xi 



D H l{ NLO) f *1 Q 2\ + D H l{ NLO) f £l q2 
q \%i J H \X\ 



D H 2 (NLO)(^ Q 2 

y Vx 2 



&s AM), , 



D H 2{ NLO) ^ Q 2^ + D H 2 (NLO) ^ q2 ^ D H l{ NLO) ^ Q 2^ 1 + 



1-Z2 



Z2 1 ^22 



X 2 



x 2 



D H l{ NLO) ( * Q A D mNLO) f *2 q2 \ + jjH 2 (NLO) f * ( * Q 

y Vxi / Vx 2 / Vx 2 / Vxi 



D H l{ NLO) t £1 ^ + ^(ATZO) / £1 Q \ -qH 2 (NLO) f *2 Q 2 
q \Xi J H \X\ J y \x 2 



x 2 t— 4 q L 2vr J Vx x 



+ DL>f H2{Ari0) f^,^,g 2 )l +2%i M )(a;) J DD^( 7Vi0 )f^,^,g 2 H. (28) 
Vxa; / 27r y y Vxa;/ 



The longitudinal contribution is obtained by replacing /( M ) — > and omitting the corre- 
sponding LO terms proportional to 5(1 — xi)8(l — x 2 ) and 5(1 — x). 

While the dependence on the momentum fractions of the double-fragmentation functions can 
not be computed within perturbative QCD, the factorization scale dependence is driven by the 
evolution equations. As for the ordinary fragmentation functions, these equations can be obtained 
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by requiring that the bare functions DDf lH2 do not depend on the factorization scale 

' / DD? lH2 (x,y) = 0. (29) 



dlogM 2 

Replacing DDj IlH ' 2 (x J y) from Eq. (j27j) we obtain: 



d 



b iH2{ nlo ){ 2) = «. /* * p (tt)fl ^p y ^ (30) 
a log M 2 27r it J \u u J 

+ 27T4 «(!-«) U * VI- 



The first term in the right hand-side corresponds to the usual homogeneous evolution of the 
fragmentation functions Df. It indicates that the probability of obtaining the hadrons Hi and 
H 2 from the parton % is affected by the possibility of the emission of a parton j with momentum 
fraction u, which can produce two hadrons by a double fragmentation. The second term, on the 
other hand, is inhomogeneous and it is not present in the evolution equations of the Df functions. 
It corresponds to the case of a the parton i that evolves emitting the partons j and k with 
fractions u and 1 — u respectively, which can also give origin to two hadrons, but now by means 
of the mechanism of single fragmentation form each one of them. Both terms are represented 
symbolically in the FigEJ These equations fully agree with the ones originally proposed in [HE). 



dDD, 
clLogM 2 




Hi / 
D J // H 2 



Figure 4: Graphic representation for the evolution of the DDf lH2 . The first diagram corresponds 
to the first term in Eg. The last one to the inhomogeneous term. 

The presence of these two terms in the evolution equations evidences the fact that, within the 
precision of any possible detector, it is physically impossible to determine which mechanism, either 
single or double fragmentation, has been responsible for the production of two hadrons when they 
are found in the same jet. In the same way, and beyond LO accuracy, only the sum of the two 
terms in Eq.© associated to each one of the mechanisms that contributes to the cross section, 
has a physical meaning. In this sense, the similarity with the situation of the fracture functions^ 
in DIS appears again. 
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4 Conclusion. 



In this work, the cross section for the production of two hadrons in e + e~ annihilation is calculated 
to order a s considering events that include the possibility that both hadrons appear in the same 
jet. For this purpose it is necessary to extend the fragmentation model including a new type of 
functions, the double- fragmentation functions DD^ lH2 , that describe the transition of a parton 
into two hadrons. These functions, along with the single-fragmentation function Df, allow an 
unified treatment for the description of two-hadron production in e + e~ annihilation. 

While at leading order the DDf lH2 functions are necessary to contemplate the possibility of 
the double fragmentation, at next-to-leading order and beyond, they are required to perform the 
factorization of divergences that cannot be absorbed in the single-fragmentation functions. As 
a result, they obey the inhomogeneous evolution equations in Eq.()30p. where the two mentioned 
mechanisms of fragmentation are involved. We showed, for the first time, that introducing the 
double-fragmentation functions the usual factorization procedure can be enlarged consistently for 
the production of two hadrons to order a s , reobtaining the evolution equations originally proposed 
in |11 El- 
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A Appendix 



The NLO(MS) corrections to the single fragmentation term are given by 



fgqB ( X l, X 2) 



+ 



+ 



3 (1 -xi)(l -ar 2 )+[o,i] 
4 



+ 



r i 4 

P OT (a;i) log |^(1 - xijxij + -(1 - xi) 



6(1 -x 2 ), 



3 (i - ^) + :L + (r- ^) +[0; . + - - - )( - 2 - s) + 



3(1+^) 



log(l - an) 
1 — #i 



/ log(l ~X 2 ) \ 

V i-x 2 ; 



+[0,1] 



+ -(1 - Xi) + P M (xi) logxi 



+[0,1] 



+ o( X - ^2) + ^(^2) logx 2 



6(1 -x 2 ) 
6(1 -x ± ), 



x\ + (2 - xi + x 2 ) 2 



3(1- xi) (xi + x 2 - l)+ [i-x 2 , 1] 



+ 



Pgg(Xl) !Og(l - Xl)Xl 



+ -Xl) 



5(xi + x 2 - 1), 

4 x 2 + (2 - xi - x 2 ) 2 

3 (1 - xi)+[ ,i] (xi + x 2 - 1) 

4 x 2 + (2 - X! - x 2 ) 2 



3 (1 - x 2 )+[ ,i] (xi + x 2 - l) +[ i- a;2 ,ii 

4 



r 1 4 

P OT (x 2 ) pog(l - x 2 )x 2 J + -x 2 

Pgq(x 1 ) log(l - Xi)xi 



+ 



5(1 -Xi), 

4 

+ 3X1 5(1 



P OT (xi) log(l - Xi) 2 Xi 

4 x 2 + (2-xi-x 2 ) 2 
3 (1 - x 2 )+[o,i] (xi + x 2 - 1) 

'x x + x 2 - 1 



-Xi 

3 



+ 



<5(xi + x 2 - 1), 







4 




log(l - Xi)xi 


+ 3*! 



5(1 -x 2 ), 



/ Xi +x 2 - 1 \ 

V Xo / ' 



fq9A/B( Xl ' X2 ) — 3 ( x 2 )' 

/ota/b ( Xl ' 3:2 ) = 3 ( ^2 ^ • 



12 



The corresponding corrections to double fragmentation are 




logx - - 



3 1 3 5 



2 (1 - aO+[o,i]+ 2 2 



/j M) (*) 



+ 



P OT (x)[log(l-a;)x 2 ] 





4 
3' 
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